Time averaged total force on a dipolar sphere in an electromagnetic field 
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We establish the time averaged total force on a subwavelength sized particle in a time harmonic 
varying field. Our analysis is not restrictive about the spatial dependence of the incident field. We 
discuss the addition of the radiative reaction term in the polarizability in order to correctly deal 
with the scattering force. As a consequence and illustration, we assess the degree of accuracy of 
several polarizability models previously established. 



In the last years there has been an increase of interest 
in the manipulation of small particles by means of the 
light Lorentz's force. For spheres subwavelength radius 
the total force due to a light wave is usually split into 
two parts from the use of the dipole approximation (cf. 
Ref. [1]): a gradient force (p.V)E, essentially due to the 
particle induced dipole moment p interacting with the 
electric field E; and a scattering and absorbing forces 
■ip x B, where B is the magnetic vector, p = dp/dt, and 
c is the speed of light in vacuum. It has been customary 
after Ref. [1] to express the gradient force F gra( j as (see 
e.g. Ref. [2]): 



F grad = (l/2)a VE 2 



(1) 



where ao is the particle polarizability, satisfying the 
Clausius-Mossotti equation: 
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a being the particle radius and e denoting the dielec- 
tric permittivity. On the other hand, the absorbing 
and scattering forces are written in the approximation 
of small spheres through the absorbing (C a b s ) and scat- 
tering (C sca t) cross sections as: 
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where k represents the light vector [k — |k|). On using 
the expression of these cross sections in the dipole ap- 
proximation, just the first term of their Taylor expansion 
versus the size parameter, x — 27ra/A, is usually consid- 
ered. 3 

At optical frequencies involved in many experiments, 
however, only the time average of the electromagnetic 
force is observed. In this letter, we establish the form 
of the time averaged total force on a particle with- 
out restriction on the spatial dependence of the elec- 
tromagnetic field. Further, we discuss some of its con- 
sequences. For time harmonic electromagnetic waves, 4 
we write E(r,t) = E e- Mt , B(r,i) = B e~ lujt and 
p(r,i) = poe~ lwt , Eo, Bo, and po being complex func- 
tions of position in space, and 5fte denoting the real part. 
Then, the time average of the total force is: 



<F> =3?l T/2 
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(p + p*).V(E + E*) 

dt, (4) 



-(p + p*) x (B + B*) 

c 



where * stands for the complex conjugate. On perform- 
ing the integral and using E , B , and p , Eq. (4) yields 
for each i th Cartesian component of the averaged total 
force: 



< F l >= (l/2)3?e 



(5) 



for (i — 1, 2, 3), where Sijk is the Levi-Civita tensor. On 
using the relations Bo = ^|jV x Eo, po = aEo, and 
Po = — iujpo one gets for Eq. (5): 



< F i >= (l/2)SRe 



(6) 



On taking into account that 
can finally express < F l > as 



-8 l m Sj one 



< F l >= (l/2)SRe [aE 0j d i (E j y 



(7) 



Eq. (7) is the main result of this letter. It represents the 
total averaged force exerted by an arbitrary time har- 
monic electromagnetic field on a small particle. 

In this connection, Ref. [5] establishes the average force 
on an object represented by a set of dipoles when the elec- 
tromagnetic field is a plane wave. We notice that in this 
case Eq. (7) reduces to just Eq. (3), in agreement with 
Ref. [5]. However, as we shall illustrate next, Eq. (7) 
permits to apply the couple dipole method (CDM) to 
more complex configurations like that of a small particle 
in front of a dielectric surface, under arbitrary illumina- 
tion (see Ref. [6] for a discussion on the CDM for large 
particles). Also, the absence of the magnetic field Bo in 
Eq. (7) eases the computations. 

Conversely, when Eq. (2) for the polarizability is intro- 
duced into Eq. (7), one obtains for the i th component of 
the time averaged optical force: 
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<F t > = (l/2)a 5Re ^^o)* 

= (l/4)a 5Re [d 4 |E | 2 ] = (l/4)a (c1 4 |E | 2 ) (8) 

which is just the gradient force. Notice the factor (1/4) 
(see e.g. Ref. [7]) instead of which the factor (1/2) for 
non-averaged fields often appears in the literature (see for 
example Refs. [2,8,9]). In agreement with the remarks of 
Ref. [10], now the scattering force, Eq. (3), vanishes and 
thus, < F > reduces to the gradient force. Therefore, 
ao must be replaced from its static expression (2) by the 
addition of a damping term. This was done by Draine, 10 
who with the help of the optical theorem, obtained: 
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ao/(l - (2/3)ik 3 a ) 



(9) 



The existence of the imaginary term for a in Eq. (9) is 
essential to derive the correct value for the averaged total 
force due to a time varying field. 

As an illustration, let the field that illuminates the 
particle be the beam whose electric vector is: 

E x = cxp(-x 2 /2) exp(i(fcz - urt)), E y = 0, E z = (10) 

On using Eqs. (2) and (10) in Eq. (7), wc find: 

< F x > = -(a /2)xcxp(-x 2 ) (11a) 

< F z > = 0. (lib) 

On the other hand, if the correct polarizability, Eq. (9), 
is introduced with Eq. (10) into Eq. (7), the total force 
is then expressed as: 

<F X > = (l/2)3?e [-axexp(-x 2 )] 
— (cno/2)a;exp(— x 2 ) 



1 + (4/9)fc 6 a 2 
<F Z > = (l/2)fcexp(-a; 2 )5Re 
exp(-2; 2 )fc 4 a 2 /3 
= l + (4/9)fc 6 a 2 ' 
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(12a) 



(12b) 



For a particle with a radius a <C A, e.g. a = lOnm, 
at wavelength A = 632. 8nm and e = 2.25, the factor 
(1 + (4/9)fc 6 a{)) is very close to one (notice in passing 
that the expression used in Ref. [11] for a makes this fac- 
tor unity). We thus see that in contrast to Eqs. (11), the 
correct form for the polarizability, Eq. (9), leads to a total 
force given by Eqs. (12a) and (12b), which can be associ- 
ated to the gradient and scattering components, namely, 
to the time average of Eq. (1) and Eq. (3) with C a bs = 0, 
respectively. 

In the case of an absorbing sphere, the dielectric con- 
stant becomes complex and so is ao- Then, Eqs. (12) 
with a <C A become: 



<F X > = -{l/2)$te[a ]xexp{-x 2 ) 
exp(-x 2 )fc 4 |a | 2 



<F Z > 



kexp(-x 2 ) 
+ ~ 3m[a ]. 



(13a) 



(13b) 



The imaginary part of ao does not contribute to the com- 
ponent < F x >, namely, to the gradient force Eq. (13a). 
On the other hand the absorbing and scattering force, 
Eq. (13b), exactly coincides with the expression obtained 
from Eq. (3). 
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FIG. 1. a) Relative difference between the force computed 
by the exact Mie calculation and by the dipole approxima- 
tion: CM-RR (full line), LAK (thick line), DB (dashed line). 
The sphere is of glass (e = 2.25), illuminated by an incident 
propagating plane wave (A =600nm). b) Same as Fig. la for 
a silver sphere (A =400nm, e = — 4 + i0.7). 
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FIG. 2. a) Relative difference between the component of 
the force perpendicular to the incident wave vector obtained 
from CDM and by the dipole approximation: CM-RR (full 
line), LAK (thick line), DB (dashed line). The sphere is 
of glass (e = 2.25) illuminated by an incident evanescent 
wave (A =600nm). b) Same as Fig. 2a for a silver sphere 
(A =400nm, e = -4 + i0.7). 

We next illustrate the above arguments with some nu- 
merical calculations that permit us to assess the degree 
of accuracy of several polarizability models previously 
established. We first compare the relative difference be- 
tween the force obtained from the exact Mie calculation 
and the most usual polarizabilities models, namely, those 
of Lakhtakia 12 (LAK), Dungey and Bohren 13 (DB), and 
the Clausius-Mossotti relation with the radiative reaction 
term (CM-RR) 10 previously discussed, versus the radius 
a of a sphere illuminated by a plane propagating wave in 
free space (Fig. la and lb). Secondly, when this sphere 
is illuminated by an evanescent wave created by total in- 
ternal reflection on a dielectric surface, the component 
of the force perpendicular to the incident wave vector 
(Figs. 2a and 2b) is compared with the result derived 
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from the CDM (as discussed in Ref. [6]). All curves are 
represented up to a = A/10. The relative difference (%) 
plotted is defined as: 100 x (F re f — F po i,)/F re f where 
pol denotes the force obtained from the corresponding 
method used for the polarizability (among LAK, DB, 
CM-RR) and ref stands for the force derived from the 
Mie calculation when the incident wave is propagating, 
and from the CDM when the incident wave is evanescent. 

We first consider a dielectric sphere (glass, e = 2.25) 
illuminated at A = 600nra (Figs, la and 2a). We ob- 
serve that, for an incident propagating wave (Fig. la), 
the result from the CM-RR relation is better than that 
of DB, and this, in turn, is better than the result from 
LAK. The force over a dielectric particle given by the ex- 
act Mie calculation is: F = C scat (l — cos 9)\E\ 2 /(8n), 
and that obtained from the dipole approximation is: 
F = (l/2)\E\ 2 $te[-ia]. When the DB model is used, 
then a = (3/2)zai/fc 3 where a\ is the first Mie coeffi- 
cient, hence, 47r5fte[— ia] is the scattering cross section 
for an electric dipole. However, when Eq. (9) for the 
CM-RR is employed, 47r5fte[— ia] constitutes only the first 
term of the Taylor expansion of the scattering cross sec- 
tion versus the size parameter x. This is why C sca t is 
underestimated when it is calculated from the CM-RR 
model. Therefore, the DB model should be better. How- 
ever, in both cases the factor cos# has not been taken 
into account in the dipole approximation and, thus, both 
results overestimate the force. Hence, this factor cos 6* 
makes a balance making the CM-RR result closer to the 
Mie's solution. In the case of an incident evanescent wave 
(Fig. 2a), DB and CM-RR results are very close together, 
this is due to the fact that the real parts of both polar- 
izabilities are very close to each other. We see that LAK 
result, as with a propagating wave, is far from the correct 
solution. 

As a second example, we consider a metallic sphere (sil- 
ver) illuminated at A = 400nm (e = —4 + i0.7). We now 
observe that for an incident propagating wave (Fig. lb), 
the DB model yields the best result. The force can ex- 
actly be written as F = (C ex t + C sca tCOs9)\E\ 2 /(8tt). We 
notice that now C sca tCOS 9 is of sixth order in x in compar- 
ison with C ex t- Since C ex t oc 3Je[ai] in the electric dipole 
limit, the DB formulation appears as the best here. Also, 
for incident evanescent waves (Fig. 2b), DB gives the 
most accurate solution. However, for a metallic sphere, 
the relative permittivity much depends on the wavelength 
used, hence, it is now difficult to establish a generaliza- 
tion of these results. We have checked, notwithstanding, 
that for a gold or silver sphere in free space in the visible, 



DB is often the best. 

In summary, we have established the average total 
force on a little particle in a time harmonic varying field 
of arbitrary form, and thus clarify its use in the inter- 
pretation of experiments, as well as in some previous 
theoretical works. For instance, we see that Eq. (7) is 
not just the gradient force as stated in some previous 
work (see e.g. Ref. [14]). Also, this general expression 
shows the importance of the radiative reaction term in 
the polarizability of the sphere put forward by other au- 
thors. Its derivation makes no assumptions about the 
surrounding environment. It is just necessary to know 
both the electric field and its derivative at the position 
of the sphere, and thus it allows an easy handling of il- 
luminating evanescent fields. An immediate important 
consequence is that it permits to assess the adequacy of 
several polarizability models. 
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